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We investigate the recent suggestions by Barausse et al. (astro-ph/0501152) and Kolb et al. 

(hep-th/0503117) that the acceleration of the universe could be explained by large superhorizon 
fluctuations generated by inflation. We show that no acceleration can be produced by this mech¬ 
anism. We begin by showing how the application of Raychaudhuri equation to inhomogeneous 
cosmologies results in several “no go” theorems for accelerated expansion. Next we derive an exact 
solution for a specihc case of initial perturbations, for which application of the Kolb et al. ex¬ 
pressions leads to an acceleration, while the exact solution reveals that no acceleration is present. 

We show that the discrepancy can be traced to higher order terms that were dropped in the Kolb 
et al. analysis. We proceed with the analysis of initial value formulation of general relativity to 
argue that causality severely limits what observable effects can be derived from superhorizon per¬ 
turbations. By constructing a Riemann normal coordinate system on initial slice we show that no 
infrared divergence terms arise in this coordinate system. Thus any divergences found previously 
can be eliminated by a local rescaling of coordinates and are unobservable. We perform an explicit 
analysis of the variance of the deceleration parameter for the case of single held inflation using usual 
coordinates and show that the infrared divergent terms found by Barausse et al. and Kolb et al. 
cancel against several additional terms not considered in their analysis. Finally, we argue that in¬ 
troducing isocurvature perturbations does not alter our conclusion that the accelerating expansion 
of the universe cannot be explained by superhorizon modes. 

PACS numbers: 98.80.Jk, 98.80.Cq 


I. INTRODUCTION 

There are now several lines of evidence pointing toward 
an acclerating expansion of the universe. These include 
the luminosity distance-redshift relation measured from 
Type la supernovae (SN la) ; the com¬ 

bination of the angular diameter distance to the surface 
of last scattering and the physical matter density 
measured from the cosmic microwave background (CMB) 
with the low values of UlmMo favored by large-scale struc¬ 
ture data 13) IE llE HE : and, most recently, the inte¬ 
grated Sachs-Wolfe effect ■ 

It is well-known that such an accelerating expansion is 
impossible if one makes the following three assumptions: 

1 ) the strong energy condition (SEC) holds, i.e. the 
density and isotropic part of the pressure seen by all 
observers on timelike trajectories satisfy /o+3p > 0; 

2 ) the universe is described by general relativity 
(GR); and 

3 ) the universe is homogeneous and isotropic, in par¬ 
ticular the Friedmann-Robertson-Walker (FRW) 
metric is applicable. 

Any explanation for the acceleration of the universe must 
drop at least one of these three assumptions. Usually ei¬ 
ther assumption m or m is dropped. In such a case, 
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we use the term “dark energy” to describe any SEC- 
violating matter field, and “modified gravity” to denote 
any departure from GR. These explanations for the ac¬ 
celeration could be considered unsatisfying since there is 
presently no other experimental motivation for modifica¬ 
tions to GR, and the matter helds normally considered 
in cosmology, including baryonic matter, photons, neu¬ 
trinos, and cold dark matter (GDM) all obey the SEG. 
In either case, new physics must be invoked. In contrast, 
it is observed that assumption dsi is not exactly valid in 
the real universe. Therefore several recent papers 
have asked whether in fact the “backreaction” from these 
perturbations to the universe can explain the accleration, 
without dropping the SEC or GR. 

The purpose of this paper is to examine the recent 
suggestions by Barausse et al. (hereafter BMR) and 
Kolb et al. (hereafter KMNR) that perturbations 
on scales larger than the Hubble length can explain the 
acceleration. In particular, these papers suggest that the 
time evolution of these large-scale perturbations produce 
a large variance of the deceleration parameter q. Since 
potential perturbations at the horizon scale are of order 
~ 10“®, one would expect the fluctuations in q to be of 
this order, however KMNR argues that corrections due 
to very large-scale modes (hundreds of e-folds outside 
the horizon) can cause the standard deviation of q to be 
^ 10“®. In particular, for spectral index rig < 1 they 
claim that the corrections from very large-scale modes 
contain an infrared divergence. If the variance is very 
large, this could mimick dark energy and cause an appar¬ 
ent acceleration. Indeed, the existence of perturbations 
on scales well beyond the horizon is likely in the context 
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of inflation [Hilil, since the inflaton fleld may have 
been in the “slow roll” regime for many more e-folds than 
the minimum of ~ 60 necessary to solve the horizon prob¬ 
lem. This paper focuses on the possibility of accelerated 
expansion from super-horizon perturbations; we do not 
consider mechanisms based on sub-horizon physics here. 

We argue that in fact no acceleration is produced. Our 
arguments can be summarized as follows: 

• The result that GR-I-SEC forbids accelerating ex¬ 
pansion is more general than the FRW cosmol¬ 
ogy. The application of the Raychaudhuri equation 
to inhomogeneous cosmologies immediately yields 
several “no go” theorems for accelerating expansion 
with zero vorticity, depending on how acceleration 
is defined. One of these theorems rules out the 
KMNR mechanism. We treat the definition of the 
deceleration parameter q and the conditions under 
which (7 < 0 is possible in Sec. im 

• In Sec. EH we examine the physical reason why 
KMNR found that their deceleration parameter be¬ 
comes negative, in contradiction to our conclusion 
in Sec. es In particular, we examine one configu¬ 
ration of initial perturbations that can be solved 
exactly; we show that in this case, the discrep¬ 
ancy arises because of higher-order terms in the 
perturbative expansion of fit (their Eq. 7) that were 
dropped in their analysis. Once the full result for 
fit is used, the deceleration parameter q never be¬ 
comes negative. The results of Sec. lIIIl are less gen¬ 
eral than those of Sec. CH since they only apply to 
a specific configuration. Nevertheless they provide 
important insight into the calculation of KMNR. 

• Causality severely limits what observable effects 
can be derived from super horizon perturbations. 
If one constructs an initial 3-dimensional hyper¬ 
surface E at the end of inflation, then any ob¬ 
servable quantities can in principle be determined 
purely from the values of the initial perturbations 
on E n J~{0), where J~{0) is the causal past of 
the observer. There may be perturbations on very 
large scales, e.g. due to the early stages of inflation, 
but one can only observe these to the extent that 
they affect the initial conditions on ’SC] J~ (O). We 
show in the first half of Sec. IIVI that the observable 
effect of superhorizon perturbations found by BMR 
- including the infrared divergence - is in fact due 
to a local rescaling of the coordinates caused by 
the particular choice of gauge for the superhorizon 
modes. 

• In Sec. we calculate the contribution of very 

long-wavelength modes to the deceleration parame¬ 
ter in single-field inflation models. We find that this 
contribution vanishes when all effects, including 
non-Gaussianity and loop corrections to the scalar 
power spectrum, are taken into account. This 


means that the long-wavelength modes in single¬ 
field inflationary models cannot produce a large 
variance of the deceleration parameter. 

Sec. El shows that regardless of the potential perturba¬ 
tions, “acceleration” as defined by KMNR cannot be 
obtained without sacrificing the SEC or GR or invok¬ 
ing vorticity (which is not present in KMNR and is 
not predicted by inflation), while Sec. argues that 
there is no infrared divergence and that the fluctuations 
in q at the horizon scale are indeed of order ~ 10“®. 
Secs. EH and iTvn deal with the more specific issues of 
why KMNR found acceleration and infrared divergence, 
respectively, and address technical points about their cal¬ 
culation. Note that our arguments involving the “no go” 
theorems are logically independent from the arguments 
involving causality. Either one rules out the KMNR 
mechanism as the explanation for the accelerating ex¬ 
pansion of the universe. 

The results of Sec. ITvTI are specific to single-field in¬ 
flation, and do not apply to other models. For example, 
in multi-field inflation large isocurvature perturbations 
are possible, which can produce large variances in the 
deceleration parameter when one averages over different 
portions of the universe |^. However the analysis of 
Sec. El still applies to these cases, i.e. for several of the 
definitions of the deceleration parameter there is no pos¬ 
sibility of obtaining g < 0 without rejecting the SEC or 
GR. 

The notation used in this paper is as follows: the 
metric signature will be —h H—h, with hij denoting the 
spatial metric. Greek indices will range over {0,1,2,3} 
whereas Latin indices range only over the spatial direc¬ 
tions {1,2,3}. In cases where we do explicit calcula¬ 
tions with the metric, we will use the synchronous CDM- 
comoving gauge, in which <700 = —1, goi = Oj and the 
CDM particles move along curves of constant spatial co¬ 
ordinates (we show in Sec.Elthat this gauge exists in 
any scenario, including KMNR, where the CDM vorticity 
vanishes). Symmetrization of indices will be defined by 
A(^ab) = ^{Aab + Ata), and similarly for antisymmetriza- 
tion, which is denoted with square brackets A^g_^. In 
the model of KMNR, there are only scalar perturbations 
present initially, and the spatial metric becomes 

at early times when the perturbation mode of interest is 
outside the horizon (but after the end of inflation); we 
use this as the definition of ip. The overdot will denote 
derivatives with respect to the physical time coordinate 
t] we will not do any calculations with the conformal 
time in this paper. The analyses of Sec. EH and Sec. Cvl 
are done on a flat FRW background with a scale factor 
a{t) oc normalized to a(to) = 1 today. 
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II. DECELERATION PARAMETER 

In order to make a statement about the universe ac¬ 
celerating, it is necessary to have a precise definition of 
the deceleration parameter that makes sense in non-FRW 
cosmologies. In the recent backreaction papers, several 
definitions have been proposed. We consider these here. 


A. Definitions 

The simplest definition is based on the congruence of 
timelike geodesics formed by the CDM particles with 4- 
velocity u^; one may define a local Hubble constant 

Hi = ( 2 ) 


and this is recognized as the definition of the variable 
a;*)” used by Ref. It is also equivalent to the 
operational definition used by KMNR since our det h is 
equivalent to their and they ignore the small- 

scale perturbation modes 'kg. In the latter case where a 
deceleration parameter is used, it is defined as 


ad^a _ dtHi _ 

(W 


(7) 


i.e. the deceleration parameter q used by KMNR is equiv¬ 
alent to our < 71 . 

The Hubble and deceleration parameters Hi and qi as 
defined above are functions of both the time coordinate t 
and the spatial coordinates a;*. It is sometimes proposed 
to measure spatially averaged values of the Hubble con¬ 
stant, 


and then a deceleration parameter 

q^ = -l-H{\>^V^Hi. (3) 


H2{t) 


fy Hi'/hd^x 
Jy '/h (Px 



( 8 ) 


Here p, is the Lagrangian proper time derivative 
d/dtproper associated with a CDM particle. Also Eq. 0 
implies that Hi is related to the expansion tensor via 
9 = 9^ = 3Hi. Since the CDM particles follow geodesics, 
the expansion tensor is purely spatial in the CDM parti¬ 
cle’s frame, i.e. it is orthogonal to u: 9fj_^u'' = 9f^^u^ = 0 
(c.f. Eq 9.2.2 of H^). We will need to decompose 9^^, 
into the usual irreducible tensors 

d/ii/ = (4) 

where g^i, + is the spatial metric (where “spatial” 
is defined by an observer sitting on a CDM particle), 
is a 3 X 3 purely spatial antisymmet¬ 
ric tensor with 3 independent components, and is 
a traceless-symmetric spatial tensor with 5 independent 
components. 

In the special case where the vorticity 2-form uj = ^du 
is zero, we can integrate the exact 1 -form u and write 
u = —dt for some locally defined 0-form t. (In most 
cases of interest, the spacetime is simply connected and t 
is defined globally.) The congruence of CDM particles is 
then said to be irrotational. If one labels the CDM par¬ 
ticles with coordinates x*, then u = —dt implies that the 
covariant components of u are Up = (—1,0,0, 0). How¬ 
ever since the particle spatial coordinates are constant, 
M* = 0 ; normalization UpU^ = — 1 requires = 1 and 
= (1,0, 0,0). The requirement that these components 
of Up and be consistent forces goo = — 1 and goi = 0 , 
and the line element of the metric then becomes 

ds^ = —dP + hij{x^,t) dx^dxP (5) 

A straightforward computation shows that 

Hi = = ^h^^dthij = iStlndeth, ( 6 ) 


and 


92(i) = 


Jy qiVhd^x 

fy 'Jh d^x 


(9) 


Note that, unlike the case of qi, we have <72 yf — 1 — 
H 2 IHI 

The deceleration parameters <71 and 52 are convenient 
because they are easy to compute in many cosmological 
models, avoiding the complicated process of tracing null 
geodesics and estimating luminosity distances used by 
BMR. The latter, of course, is their principal shortcom¬ 
ing: qi and (72 do not correspond to the observables in 
any of today’s major cosmological probes. The famous 
Hubble plot of z versus dz, allows one to measure Hi, 
however extragalactic astronomy is only ~ 10 ^ years old 
and this is a woefully inadequate baseline with which to 
directly measure the time derivative Hi and hence the 
deceleration parameter qi. We also cannot cover the en¬ 
tire universe with astronomers so as to measure spatially 
averaged quantities such as H 2 or 52 - 

Since SN la luminosity distances are one of the major 
pieces of evidence for the accelerating expansion of the 
universe, one way to make contact with observations is to 
calculate the luminosity distance-redshift relation. This 
was first done in the context of backreaction by BMR, 
who expanded the luminosity distance dz as a Taylor 
series in redshift. 


dL = H-,^z+^^^z^ + 0{z^). ( 10 ) 

As defined here, the Hubble and deceleration parameters 
ids' and qo' depend not only on the observer’s spacetime 
coordinates t, x*, but also on the direction of observation 
(RA and Dec), because in a perturbed universe there is 
no reason why the dL-z relation should be the same in 
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every direction. BMR then went on to average these over 


solid angle to obtain 


H 3 = {H3>)at: and 93 = {q3')A7v 

( 11 ) 

A related procedure is to write Eq. 11 OB as a 
z{dL), 

function 

z = H3'dL — -773/(1 — (73')^i + 0{d\). 

( 12 ) 

This motivates the definition 


{z)a. = HAdL - ^ 4(1 - 94)^2 ^ 

(13) 


where we take the angular average over the redshift at 
fixed c?L rather than averaging the Taylor coefficients. 
This procedure has the advantage of making sense in 
slightly more general conditions. In particular, if dy 
has a direction n*, not necessarily an eigenvalue, where 
Oijn'"n^ = 0 , then Hy = 0 in that direction and 53 / in 
Eq. (IIOII is undefined. Then <73 is not defined either. In 
contrast, Eq. OS breaks down only under the more re¬ 
stricted condition that H 4 = 0, i.e. that the angular 
average of vanishes. As we will see in the next sec¬ 
tion, it is also easier to prove theorems about <74 because 
one can avoid angular averages with quantities such as 
Hy in the denominator. 

One can relate and (74 to and <73 by considering 
the angular average of the coefficients oi and d\ in 

Eq. itT^ . 

Ha = {H 3 ')a 7 , = H 3 and 
Hia-qA) = (i/f,(l-<73'))4.. (14) 

Thus we always have H 3 = Ha- If H 31 is isotropic, then 
it cancels out in the second line implying qs = Qa, but 53 
and qA need not be identical in general. 

B. Conditions for acceleration 

Having written down several definitions of the deceler¬ 
ation parameter, one can ask: under what conditions can 
the universe be accelerating, i.e. under what conditions 
is (7 < 0 possible? 

The condition for acceleration is easiest to derive in the 
case of < 71 . We begin with the Raychaudhuri equation, 

dO 0^ 

-^ = ^ (15) 

where is the Ricci tensor. Using this Eq. © can be 
re-written as 

HUi = (16) 

The Einstein equation allows us to re-write the last term 
in terms of the matter content in the CDM particles’ rest 
frame, 

= ^(o-^cr;; - -k ^(p-k3p), (17) 


where p = is the total density and p = is the 
isotropic pressure. Since a is purely spatial in the local 
Lorentz frame of a CDM particle, its square is semipos¬ 
itive definite, cr^cr)) > 0 . Thus we immediately see that 
if the strong energy condition (SEC) is satisfied (so that 
p + 3 p > 0) and there is no vorticity (w = 0), it follows 
that gi > 0. 

The vorticity uj plays a key role in all of the “no go” 
theorems regarding accelerated expansion, as well as in 
the existence of the synchronous CDM-comoving gauge 
(Eq.|SJ) and the spatially averaged deceleration (72, so it is 
worth exploring when it can be nonzero. It can be shown 
(see e.g. Eq. 9.2.14 of |^) that if a; = 0 at one time 
then it remains zero at all future times (at least until 
orbit crossings appear). In particular, the existence of 
the synchronous CDM-comoving gauge implies u = —dt 
and u! = du = —d^t = 0, so if the synchronous CDM- 
comoving gauge exists at early times then oj = 0. KMNR 
explicitly use the synchronous CDM-comoving gauge in 
their calculations, so their model has zero vorticity. This 
means that the “no go” theorems that <71 > 0 applies 
to KMNR; since they use <71 as their definition for the 
deceleration parameter, the KMNR mechanism cannot 
produce acceleration. A corollary is that <71 — > — 1 in 
KMNR must be an artifact of the perturbative expansion. 
We will return to this point in Sec. ITTTI 

Since <72 is a spatial average of qi, it follows trivially 
that in the absence of vorticity or SEC violation, <72 > 0 . 

The deceleration parameters <73 and (74 in terms of lu¬ 
minosity distances can be related to observations. The 
simplest way to do this is to consider the photon trajec¬ 
tories that arrive at the observer O, and find the lumi¬ 
nosity distance dL{v) and redshift z{v) as a function of 
affine parameter v. The photon is assumed to arrive at 
the observer at affine parameter u = 0 with unit energy, 
i.e. its 4 -momentum k satisfies = —1 at u = 0 . 

We will assume that the voriticity of the CDM parti¬ 
cles vanishes so that we can use the synchronous CDM- 
comoving gauge (Eq.[ 3 ). In this gauge, the “ 0 ” direction 
coincides with the CDM particle 4 -velocity, so the purely 
spatial nature of the expansion tensor can be used: 
^00 = Ooi = Oio = 0, and we will henceforth simply write 
Oij. We also define the photon’s 3 -velocity at the observer 
ni = ki{v = 0), and note that this is a unit vector since 
the photon is massless. Note that u < 0 on the observer’s 
past light cone. 

The redshift associated with a particular point on the 
photon’s trajectory is simply \+z{v) = fc°(u), since fc°(u) 
is the photon’s energy as measured by the CDM particle 
and the observer sees energy 1 . We can expand k^{v) to 
second order using the geodesic equation, 

dlc^ 

— = and ^ = -r“ ( 18 ) 

dv dv ^ 

The derivatives of are then 

^ = -T^kV = -ikikV 


(19) 
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since r°o = = 0 in the synchronous gauge, and 

1 

Plugging in the Christoffel symbols yields 

dkP d?k^ 


( 20 ) 


k^ = 1 


dv 

1 ■ 


du 2 2 


= 1 - hijrdrdv + 

2 

hi-jU^ 


'{hij Tl 


k^m 


+ ^^(2h*'=hfczn' - h^%mk,in^n 




( 21 ) 


Introducing the extrinsic curvature Kij = —^hij, and its 
3-dimensional covariant derivative 

k^{ij\k) ^{ij,k) d” , (^^) 

yields the redshift 

2 

z = Kijrdn^v + (^Kij + + AK^Kkj^ rdn^^. 

(23) 

Since u is the unit normal to the hypersurfaces of con¬ 
stant t, it follows that Kij = —Otj (c.f. Eq. 9.3.19 of 

[ 23 ). 

The luminosity distance to a spherical source whose 
surface is at affine parameter Vs is given by Eq. (2.14) of 
Ref. m 


dL(ys) = 


(24) 


where A is the photon amplitude (normalized by = 
A^k^k'^) and R is the radius of the source; the amplitude 
A is given by 


A oc exp ( — 


1 


9 dv 


(25) 


where 9 = is the photon expansion. The photon 

expansion near the source can be expanded as 


9 = 


V - (Va + AVs) 


+ 0[v - {Vs + Aus)], 


(26) 


where —Avs is the difference affine parameter from the 
source’s center to its surface at radius i?; this can 
be shown by the same argument leading to BMR’s 
Eq. (A27). This affine parameter can be found by not¬ 
ing that k^ = dx^/dv. Here the photon’s energy at the 
source is = 1 -I- z and the time it takes to travel from 
the center to the surface is i?, so 


, /c° 1 + ; 


(27) 


Also, integration of Eq. (1^ using Eq. (EH gives 
A(us) -Avs 1 + z 


(28) 


A(0) -Vs-Avs -VsR' 

where the last equality uses |Aus| <C |us|. This yields 

c^L('ys) =+ 2(us)]^-t-O(u^). (29) 

Comparing to Eq. (E3 gives the result, to second order, 
that 


z = — 


{l + zf 


■ + + 4.K/Kkj) 


(30) 


or 


z = —KijKn^dL — 2 {Kijrdn^Yd\ 


+ + AK/Kkj) (31) 

We can then calculate the deceleration parameter by 
comparing Eq. Ijdl |l to Eq. HI 311 . The ATyife term drops 
out in the angular averaging, and we have {rdni)^^^ = .j, 
so we get 


Ha = -Ik 


(32) 


where K = K], Note that this implies H 4 = Hi. The 
deceleration parameter is 


94 — 1 + 


-2 


/ilbTfcO 1 . 

-AK^.Kki - - -+ -K.jK^ 


+ -K"K^ 
3 * 


(33) 


This can be simplified by separating K into its isotropic 
,K/ = -H, 

ki,K^ 4. 


and anisotropic pieces, K/ = —H 46 j — cr*. Then we find 


94 = 1 


3i7? ' 5^^ 


(34) 


A further simplification comes from the formula for Kij 
(Eq. 21.82 of Ref. jH), 

R\oj = -ik,j + K,kKfy, (35) 

the antisymmetry property of the Riemann tensor cou¬ 
pled with the Einstein equation allows us to write 

K^R°aOj =R°o = SttG (t° - = -47rG(p + 3p), 

(36) 


hence 


K^jK^ = AttG{p + ip) - mi - a)al. (37) 


Substitution into Eq. El yields 


94 = 


47rG(p + ip) 7cr*(7^ 


mi 


imi ■ 


(38) 
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Once again, we see that if the SEC holds then g4 > 0 . 
Thus by this definition too a universe with only matter 
cannot accelerate. 

The deceleration parameter 53 is only well-defined for 
a positive-definite spatial expansion tensor Oij. It is also 
more difficult to work with than 54: from Eq. m , we 
find 


^ , I r, id'^z/ddl)\dL=o 
^ {dz/ddLU=or' 


( 39 ) 


The presence of variables in the denominator of this equa¬ 
tion makes it very difficult to prove theorems about the 
angular average = (q3')47r- We do know from Eq. limi 
that H4 = H2, ■ Since the mean of the square is always at 
least the square of the mean, this implies 

Hi < (i?|,) 4 .. ( 40 ) 

Combining Eq. iBUI) with the second line of Eq. la yields 
(i^|'93')4. = {HI,)4^ - Hi + Hiq4 > Hiqi. ( 41 ) 


If the SEC holds then 94 > 0 and it follows that 
(i/|,93')4^ > 0 . This inequality does allow for the pos¬ 
sibility that 53 = {qy)ATT < 0 . This can happen if there 
is an anisotropic expansion tensor (i.e. cr 7^ 0) so that 
{Hi,q3')4Tr 7^ Hi{qy)4T^ and the universe is “accelerating” 
with 53/ < 0 in some directions rf. Nevertheless even for 
an anisotropic expansion, the inequality (I?|/53')4.n- > 0 
implies that if there are accelerating directions (53' < 0) 
then there must be decelerating directions (53/ > 0) as 
well. 


III. ACCELERATION OF UNDERDENSE 
REGIONS? 

KMNR concluded that an underdense region of a post¬ 
infiationary universe filled with only dark matter can ap¬ 
pear to accelerate (according to 51) due to the influence of 
super-horizon perturbations. However, with CDM only 
the SEC should hold, and after inflation there should 
be no vorticity. In Sec. irm we showed that the SEC 
and zero vorticity imply 51 > 0 . These conclusions are 
obviously inconsistent. The purpose of this section is to 
show, in a particular case that can be solved exactly, that 
the reason for the inconsistency is that KMNR dropped 
higher-order perturbative terms in their gravitational po¬ 
tential These terms, when incorporated into their 
Eq. ( 10 ), result in a deceleration 51 that never becomes 
negative. 

The particular case that we shall consider here is the 
potential 


<^ = (1 - \c\x? 


( 42 ) 


is sufficient by causality for qy to take on the above value 
inside our horizon. Also, by Birkhoff’s theorem, if the 
potential is spherically symmetric and takes on the form 
of Eq. for |a;| less than some radius R, the future 
evolution at \x\ < R is as computed below. 

Since for C > 0 we have < 0, the potential 

Eq. ||32) can be used as a model for an underdense re¬ 
gion. Not all underdense regions look like Eq. 103), how¬ 
ever this potential can be solved and it is conceptually 
useful for understanding why the higher-derivative terms 
dropped in KMNR are in fact important. The spherically 
symmetric underdensity and overdensity have been stud¬ 
ied extensively as models of voids in large-scale structure 
1 ^ and the formation and growth of collapsed objects 
|3n| respectively, however these studies were not aimed 
at understanding backreaction. Our aim here is to un¬ 
derstand the implications of the spherical underdensity 
for the KMNR mechanism. 

Substitution of the potential in Eq. 03 into Eq. o 
yields initial conditions for the spatial metric, 

(1 ~ i^Fr) 

This is simply the metric for a 3-dimensional open uni¬ 
verse with spatial curvature —C and a coordinate sys¬ 
tem given by the stereographic projection, analytically 
continued to the case of negative curvature. The future 
evolution of this universe is equivalent to an open FRW 
universe and thus can thus immediately be written down: 

ds^ = —dt^ -\ - ^ , (44) 

( 1 -iW)' 

where aopen(0 satisfies the Eriedmann equation for an 
open universe, which has the parametric solution 

Oopen = sVC (coshry - 1) , 

t = S(sinhr7 — 77) ’ ^ 

where S is a constant and rj is the parameter. Com¬ 
parison of Eqs. 03 and 03 shows that at early times, 
ttopenit)/ a{t) —> I. Since the evolution of a{t) is simply 
that of the Einstein-de Sitter universe, 

3 

^Hotj , (46) 

we can use this as the early-time limit of aopen(i) to de¬ 
termine S, 

S = \hIC-^I'^. (47) 

KMNR defined the metric perturbation variable dt ac¬ 
cording to 



ds^ = —dt^ + a‘^{t)e ’*^dx''dx^. 


where |a;p = {x^)'^ + + {x^)"^ and C is a constant. 

In order for the conclusions derived below to be valid, it 


( 48 ) 
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Comparison of Eq. (Ell with Eq. Ell then allows us to 
determine the full nonperturbative evolution of 'h for the 
potential tp of Eq. (El : 


(49) 

The first term in this result is independent of time and, 
as correctly argued by KMNR, it does not affect qi. The 
second term can be evaluated by expressing rj in terms of 
t and hence a{t) in the second line of Eq. II45II and then 
finding Oopen using the first line. The functional form 
ri{t) cannot be expressed analytically, but a power series 
can be developed via reversion of series. 



Substitution into the Taylor expansion of aopen(i) = 
E\/C{coshr] — 1) gives 


^ 2 ,,4 


,(t) = a{t) + —pja (t) - 


5Hi 
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a (t) 




(51) 


plugging in Eq. (EH) yields 

4.(x-. t) = In (l - - ^a(() + 

(53) 

The constant in this equation agrees with Eq. (P|l . The 
first-order perturbation terms (i.e. oc C) also agree. Be¬ 
yond first order in (7, however, two discrepancies appear. 
One is the spurious term in Eq. proportional to 
that has spatial dependence. We will consider 4/ and q 
the spatial origin = 0, where the spurious term van¬ 
ishes. [This is also where the KMNR calculation is most 
likely to be valid, since they neglected terms containing 
(V(/?)^, and = 0 at the origin for our specific poten¬ 
tial, Eq. lEl .] The other discrepancy is that Eq. (El 
is missing the higher-order spatially independent terms 
proportional to (Ca)^, (Ca)^, etc. The (Ca)^ term has 
2 n derivatives in it and presumably are not included in 
Eq. (7) of KMNR because they dropped “higher deriva¬ 
tive” terms. 

We finally consider the implications of the difference 
between Eq. (El and Eq. (El for the determination of 
the deceleration parameter qi. This is given by Eq. (6) 
of KMNR, 


We are interested in the expression for Oopen as a pertur¬ 
bation series in C, since (p ^ 0 when (7^0. Eq. EU 
provides this series, although it turns out to also be a 
power series in a. [The reason for this is that C has units 
of length"^, and Hq is the only other quantity contain¬ 
ing units of length once we have eliminated po using the 
Friedmann equation. Thus C and Hq appear together in 
the combination C/Hq when expanding Uopen- Further¬ 
more, there is an invariance associated with the choice of 
epoch of the observer: if one places the observer at a later 
time instead of t, then the scale factors a and Oopen 
are multiplied by (since a oc Hq is multiplied 

by and the change in definition of comoving length 
multiplies the “curvature” constant C by hence 

C/Hq is multiplied by The invariance with respect 
to ^ forces the terms in the above expansion to take the 
form {C/H^y. Thus the expansion in C turns out 
to also be an expansion in a.] Substitution into Eq. (El 
gives 


gi = -1 + 


I -k 77-2# 
(l - £7-1#)^ 


(54) 


Here we have dropped the subscript on 'I'f since 'k = 
at the level of approximation in KMNR, and H represents 
the unperturbed Hubble constant. According to Eq. El . 
the value of qi at the spatial origin cc® = 0 is 


(KMNR) = -1 -k 2 + 10^0 Ca 
^ ’ (1 + i77o-2Ca)2 ^ ^ 

According to the full nonperturbative calculation, and 
defining Hopenit) = dopen(7)/aopen(7), we find 


# = H(t) - Hopenit) and 

# = 77(t)-77open(7). 

Substituting into Eq. (El yields 


(56) 


«'(x®,t) = In ^1 - - ^a(t) 

1 • 3(^2 

+ ^^a^{t) + 0{C^H^^a^). (52) 

Aside from the irrelevant time-independent term, 4/(a:®, t) 
depends only on the time. Technically the perturbation 
series in Eq. (El has a finite radius of convergence of 
\CH^^a\ < (|^)2/3 (see Appendix IXl. but if one had the 
full perturbative expansion one could analytically con¬ 
tinue it to later times. 

Our result for the evolution of the gravitational poten¬ 
tial should be compared to Eq. (7) of KMNR, which after 


where 77 (t) completely drops out if one uses the identity 
77 = — |772. We can use Eq. (El to find Hopenit), which 
gives 


and 


77open(7) — A 


sinhr; 
cosh 77 — 1 


<7i (true) 


1 

1 -k cosh rj 


(58) 


(59) 















Equations 113 and (E3) have some features in com¬ 
mon: they both approach 1/2 at early times, and as the 
universe expands they decrease, as expected for an under- 
dense region of the universe. However the long-time be¬ 
havior is different: as t ^ oo, gi(KMNR) ^ —1, as noted 
by KMNR, whereas rj ^ oo and hence ( 7 i(true) —> 0. We 
thus see explicitly that no acceleration is possible in this 
case and that the result found by KMNR is a consequence 
of neglecting higher order terms in their analysis: if the 
full perturbation series had been calculated (and analyt¬ 
ically continued to late times), KMNR would have found 
no acceleration. 


IV. SUPERHORIZON ADIABATIC 
PERTURBATIONS 

BMR have recently analyzed the effect of superhorizon 
adiabatic perturbations on the dL-z relation. In particu¬ 
lar they considered the value of the isotropically averaged 
deceleration parameter qs {{qo)Q in their notation). They 
found a large contribution to the variance of qs coming 
from the interaction of infrared and ultraviolet modes. 
In this section we re-derive this result by selecting an 
appropriate new coordinate system. This method illu¬ 
minates the physical origin of the large Var q^ found by 
BMR, namely that BMR assumed the ultraviolet modes 
to be statistically homogeneous, isotropic, and Gaussian 
in their coordinate system with the cosmic power spec¬ 
trum Ptp(k). Because of the infrared perturbations, this 
coordinate system differs from the locally FRW coordi¬ 
nate system erected by an observer, and the power spec¬ 
trum of ultraviolet modes seen by a specific observer dif¬ 
fers from P^{k). In Sec. II V Gl we consider the tp field 
generated by single-field inflation, and show that in this 
case the BMR infrared divergences are canceled by non- 
Gaussian features and loop corrections to the power spec¬ 
trum. 

The KMNR mechanism makes use of this infrared di¬ 
vergence. Specifically, they argue that in some models 
of inflation, there is a very large variance of the decel¬ 
eration parameter qi resulting from the interaction of 
infrared and horizon-scale modes. Within the context of 
single-field inflation, or indeed any mechanism in which 
small-scale perturbations are laid down after the large- 
scale perturbations have exited the horizon and frozen 
out, there is then no infrared divergence. We discuss the 
implications for KMNR in Sec. mni 


A. Gauge transformation; causality 

The initial state of any system in classical general rel¬ 
ativity is described by a spatial metric hij on the ini¬ 
tial hypersurface, its time derivative, i.e. extrinsic cur¬ 
vature Kij , and the appropriate variables Xa to describe 
the matter fields. These quantities are not all indepen¬ 
dent since they must satisfy the energy and momentum 


constraints. It is also well-known that two distinct sets 
of initial conditions {hij, Kij, Xa\ may actually corre¬ 
spond to the same spacetime because they differ by a 
gauge transformation, which may involve either (i) a re¬ 
parameterization of the initial hypersurface, or (ii) se¬ 
lection of a different spacelike hypersurface as the initial 
slice. We will focus our analysis on using gauge transfor¬ 
mations of type (i) since (ii) is more difficult to analyze. 

If one considers an initial hypersurface S after the end 
of inflation, one need only specify the initial data on the 
portion of E that can be seen by the observer, O, i.e. 
one needs initial data on En J~{0). In the unperturbed 
standard GDM model with the usual initial hypersurface 
(a surface of constant cosmic time), En J“(C>) is a closed 
ball with comoving radius less than the horizon size Xh, 


Xh 


da 

Jq a?H{a) 


(60) 


where we have used i?(a) oc In a non-FRW uni¬ 

verse the nature of E n J“ (Olis more complicated. Nev¬ 
ertheless, Cauchy stability implies that, for suffi¬ 
ciently small perturbations, the observables seen by O 
can be uniquely determined from the initial data on any 
open set U in Y, with D E n J“((!I,FRW), where 
J“(C>,FRW) is the causal past of the observer assuming 
the FRW metric (see Fig. ^|. In particular this means 
that given such a U, one can calculate the observables 
seen by O to any order in perturbation theory. It is also 
easy to construct U: any open ball centered at the origin 
and extending out to coordinate radius Ru > Xh will do. 
We assume Ru = 2xh for definiteness. 

The initial hypersurface considered by BMR has a spa¬ 
tial metric 


= (61) 
and extrinsic curvature 

Kij{ts,x'") = -^h,j{tY:,x’^) 

= -H{t^)h,j{t^,x'") + 0 {K,jt‘^^^),{ 62 ) 

where we write 0 {Kijt^/^) to mean that the fractional 
correction is of order Instead of a;* coordinates we 
want to construct a coordinate system which is explicitly 
locally inertial. We will therefore impose a gauge trans¬ 
formation to new coordinates y^{x^) such that at early 
times, hij is very close to a^dij (without the conformal 
factor). One realization of such a coordinate system is to 
construct geodesics in all directions from the center and 
decompose their tangent vectors in terms of a chosen ba¬ 
sis: the resulting coordinates are the Riemann normal 
coordinates. 

We begin by splitting ip in Eq. into a large-scale 
and a small-scale piece, 

(fix'") = tpL{x'') + ips{x"). 


(63) 
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observer 

FIG. 1: The regions of spacetime described in the text. Initial 
data are supplied on the Cauchy hypersurface E. Causality 
implies that the observables seen by observer O can depend 
only on the initial data within the shaded region E n J~ iO). 
Within the context of perturbation theory around an FRW 
spacetime, the perturbations to the observables can be de¬ 
termined given initial data within any open set lA containing 
E n J“(C1, FRW). Superhorizon perturbations can only be 
observed to the extent that they affect the initial data inside 
U. 


We will work to first order in Lps and in the wavenumber 
kL of the large-scale modes. That is, any terms contain¬ 
ing (p| or k\ipL ~ VV(/ 3 l will be dropped. Note that 
the split into “long” and “short” wavelengths here differs 
from the split ik = -|- dts in KMNR: we are including 

both, and are doing the separation so that k^ and (ps 
are the expansion parameters, whereas in KMNR Tg is 
dropped. 

Next we construct a Riemann normal coordinate sys¬ 
tem on the initial slice, centered at a;* = 0 , with respect 
to the large-scale metric 7 ^. The construction of such a 
coordinate system begins with the choice of an orthonor¬ 
mal basis of three vectors {(ea)*}i=i tangent to S at the 
point (ti;, 0, 0,0). (We use the barred indices a to denote 
the components associated with these basis vectors and 
with the y coordinate system that will be defined based 
on them.) Such a basis is not unique; here we choose 

(64) 

where one takes the — 1/2 power of the 3x3 matrix 
of covariant components of The Riemann normal 

coordinates are then constructed via 

a;*(yJ) = expo jefy-’l, (65) 

where exp represents the exponential geodesic mapping. 
In our case, the argument of the exponential simplifies to 
the vector 

(66) 


The geodesic equation of the 7 ^ metric can be found from 
its Christoffel symbols, 


(Px^ 


j dx^ dx^ 


6 10v>i/3[ ^V’L 

dx^ 

3 1 dx^ 

dX 


d(pL dx^ dx^ 
dx^ dX dX 


(67) 


Solving this equation to second order in y® subject to the 
initial conditions x(A = 0) = 0 and x'(A = 0) given by 
Eq. lIHHll gives 


x^ = x^{X = I) 

_ 10.^1, (0)/3i , 6 20v>_l(O)/3 

y ^ 6 


+2y'-y^(pL,j{0) 


-|2/lV.(0) 


( 68 ) 


The metric in Riemann normal coordinates defined 
with an orthonormal basis is at 0 and its first spatial 
derivatives vanish (see e.g. Sec. 11.6 of [ 23 ). Thus 

iL Tjiy’") = % + OiklpLlyl"^), (69) 

where k\ipL is the second spatial derivative of pL, which 
is suppressed for long-wavelength perturbations. Note 
that within the set hi, |yp is at most of order The 

full metric is then 

+ OiklpLly^. (70) 

We may now calculate the deceleration parameter in 
the new coordinate system. BMR has found the deceler¬ 
ation parameter to be, to first order in p, 

q3 = l + lH^^yip + 0{p^). (71) 

Since Eq. CZOl) is now simply the Euclidean metric with 
a first-order perturbation, we can apply Eq. m in the 
y coordinate system to yield 

<73 = i + 0{pl,H^^klpL). (72) 

The order of the error is now since Eq. had 

error 0 (k\{pL\y\^) and causality requires that we can only 
observe out to distances |y| < Ru ~ 0{Hq^). Thus in 
Riemann normal coordinate system, which is constructed 
locally and does not make reference to anything farther 
from the observer than the horizon, the result for the de¬ 
celeration parameter is extremely simple and large am¬ 
plitude long wavelength modes have no impact on the 
deceleration parameter. 
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B. Deceleration parameter in x coordinates 

BMR expressed the deceleration parameter 53 in the x 
coordinate system rather than y. We have already shown 
that the analysis in y coordinates results in Eq. [721 and 
that superhorizon modes cannot affect the deceleration 
parameter. Nevertheless, we wish to compare our result 
to BMR, who did the calculation in the x coordinates. 
Since the final result for any observable cannot depend on 
the choice of coordinates we must obtain the same answer 
as suggested by Eq. |Z21 We will show below that this is 
indeed the case and is a result of a precise cancellation 
between several terms, some of which were not considered 
by BMR or KMNR. 

We can recover all terms in the deceleration parameter 
up to order by using Eq. to turn Vy </55 into 

derivatives with respect to x: 

V>s = (73) 

Plugging into Eq. yields 

+gl0c,^/3v2^g ] 


where T^p is the connected trispectrum of the potential in 
the X coordinate system. BMR argues that since inflation 
produces Gaussian initial conditions, the trispectrum 
term drops out, and if rig < 1, where P,p{k) oc k"‘‘~^, 
there is an infrared divergence. We show in the next 
section that this divergence disappears in the context of 
single-field inflation if one calculates the variance of 
instead of just a single term. 

C. Infrared divergences in single-field inflation 

The computation of Var((/3VxV5) is an important step 
in determining the observable deceleration parameter. 
Nevertheless, we cannot observe directly, rather 

what we observe is 93 . Therefore we want to know 
Varg 3 . This section examines the contribution of the 
very long-wavelength modes to Var predicted by single¬ 
field inflation. We are especially interested in the case 
of rig < 1, where the variance /_i =Var(p is infrared- 
divergent: I-I —>■ -too. 

In this section, we retain the terms from Eq. ( 171 ) that 
have no factors of k^, i.e. we drop and are no 

higher than second order in which is the order at 
which BMR found an infrared divergence. BMR argued 
that the <PL,i^s,i term is subdominant and in any case the 
variance of PLp has no infrared divergence for rig < — 1 . 
The equation for (73 is then 


1 

2 1 

5 _2 

25 

1 rr-2 1 

75 

50 250 2 \ 

2 

+ h!\ 



93 - 2 + ^0 ( 

k9 + 

+ JJ-PlJ 


25 .,2 


(74) 


where in the last line we have worked to order . This 
should be compared to Eq. (36) of BMR, which reads 


93 = 


"0 


5^,2 25,^ ,2 


In order to determine the mean and variance of qs we 
need a measure. The two most obvious measures are the 
coordinate measure d^x and the measure dN associated 
with the number N of CDM particles. Since at early 
times the CDM density is constant in the synchronous 
comoving gauge, it follows that dN is proportional to 
the volume element at early times, 

dN (X (78) 


25 


0 {H-% 


(75) 


We can unify these two measures by working with the 
measure d/i defined by 


In particular, the infrared-divergent term p'S/'^p identi¬ 
fied by BMR appears in Eq. m- The PL,ifs,i term is 
not infrared-divergent (unless rig < — 1 ), nevertheless it 
is first-order in kL so it appears here; the coefficient is 
double that in BMR because of the factor of 2 in the bino¬ 
mial expansion of The variance of the divergent 

term is 


YaiipNlp) 


P^k) 


d^k 


(27r)3j 


k'^Ppik) 


d^k 


(27r)3 

^ ^ d3k- 

fc3^fciT^({k,})d(^k,)n7^, 

i=i ^ ^ 


(76) 


dp = 


l + ^vjp^-x) 


d3x = J d3x, 


(79) 


where Vj are constants. The choice Vj = 0 corresponds to 
the coordinate measure, and Vj = (—5)^/j! corresponds 
to the CDM particle number measure, up to an overall 
constant that does not affect the averages. The function 
J is the Jacobian. 

In order to compute the mean and variance of (73 in 
Eq. II77I) . we must compute some averaged values involv¬ 
ing pl and 'S/'^ps- The values linear in ps that we need 
are 


= 0 , 


(80) 
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which holds by spatial homogeneity; 

{ipL^^‘Ps)yi = 0 , 


(81) 


which results from and ips having different wavenum¬ 
bers; and 

= 0, (82) 


which results from the small-scale wavenumbers being 
very large compared to the large-scale wavenumbers so 
that the triangle inequality cannot be satisfied. Note that 
these averages are all taken with respect to the coordinate 
measure (we convert to dp, later). 

We will also need the average values (v5i(Vx‘/?s)^)x 
that are quadratic in (ps- To write these, we follow 
KMNR in defining 



where P^\k) is the power spectrum of ip calculated 
by standard first-order perturbation theory in inflation 
(i.e. with only quadratic terms in the perturbation La- 
grangian). In what follows we will consider only the long- 
wavelength contributions to /_i, as these are the contri¬ 
butions that produce the infrared divergence reported by 
BMR (and in any case they dominate for Ug < 1). The 
first-order theory power spectrum is given explicitly by 
(see e.g. Eq. (2.20) of Ref. [13) 



288Tr^GH^ 

25fc3(^2 


(84) 


where if* is the Hubble rate when the mode k becomes 
super-Hubble and 0* is the rate of change of the inflaton 
field 0 at that time. (Note that our (p is denoted by — 
in Ref. jH.) 

The easiest average value to compute is 
(<pi(Vx^s)')x = (<pi)x((Vx<^s)')x = I-lh, (85) 


which results from decoupling of the long- and short- 
wavelength modes. The integral I 3 is dominated by short 
wavelengths, and (for Ug < 1) I-i by long wavelengths. 
The average values with 0 or 1 power oi tpL are how¬ 
ever more complicated because of corrections for the non- 
Gaussianity of ip. The average value with 1 power oi ip l 
and 2 of ps is determined from the bispectrum configu¬ 
ration of p where the Fourier-space triangle has one short 
and two long sides, k^ /c 2 . This configuration has 

been computed by Ref. [13 (see Eq. 4.7 in that paper) 
and is 

(<y5ki</5k2V5k3) = 

( 86 ) 

where n^^\ki) is the scalar spectral index. |36| This im¬ 
plies 

(<Pl(Vx<P5)^)x = ll-i J [n®(fci) - l]fc4pW(fc)^. 

(87) 


The average {{'^xPs)‘^}x at hrst glance appears to be 
just I 3 . However we require the lowest-order correction 
due to long-wavelength perturbations. Thus one must re¬ 
turn to the derivation of Eq. (ISSll in order to obtain the 
correct result. The argument provided by Ref. [^ for 
this situation is that the small-scale modes ps are unaf¬ 
fected by large-scale modes except that these large-scale 
modes re-scale when ks becomes super-Hubble. In par¬ 
ticular, for fixed large-scale perturbation p^, the small- 
scale modes exit the horizon (l = —^‘Pl e-folds of in¬ 
flation earlier than in the unperturbed case. Therefore 
their power spectrum, conditioned on p^, is 




(fee 


-a 13 


27r2 ^ 

(^g5¥>i,/3)3 




27r2 


pW^ke^VL/S)^ 


( 88 ) 


If we want the power spectrum of the small-scale modes 
averaged over all regions of space, then one needs to av¬ 
erage Eq. over the probability distribution of pL- 
Taylor-expanding Eq. ilSSI) to second order in p^ and do¬ 
ing the averaging yields 






+ 


{P 


i)x 


dp 


L L 


{ke^VL/^Y 


27r2 




(89) 


where the second derivative term comes from the vari¬ 
ance of pl, and there is no first derivative term because 
{pl)x = 0. The variance of pL is {p\)x = I-i- If we 
define the variable 


l(0) 


(fc) = 


k^p!p\k) (dlnfc)' 


,[fc3p(o)(fc)] 


= a?) 


{k) + [ni°)(fc) - 1] 


(90) 


which is related to (but different from) the scalar running 
a^s\k), Eq. lIS^ simplifies to 

((VVs)") = / A:V®(fc) 

(91) 

Physically, Eq. iPUll represents the “smearing” of the re¬ 
lation between k and physical scale due to the variance 
of This smears the power spectrum Ptp(k). Thus it 
is not surprising that Eq. (inij contains a second deriva¬ 
tive of the power spectrum. The smearing effect can be 
thought of as a loop correction to the scalar power spec¬ 
trum in the sense that it involves an integration over an 
undetermined momentum (that of pL , which is packaged 
into I-I here). 

The above results are most easily expressed if we define 
the integrals 


1 25 d-^k 

(fc)/_i 


(92) 
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In particular, 


The denominator is 


h,i = I k^in^^Hk) - nPi°Hk)-^ (93) 

and 

/ 3,2 = J k^af\k)P^^\k)^. (94) 

Then the average values containing (V^tps)^ are 
((VVs)^) = /3 + 

((/3z,(VVs)^) = and 

= hl-i. (95) 

Finally we come to the issue of the mean and variance 
of < 73 . The mean of any quantity X with respect to the 
measure /r is related to the Jacobian J via 

(V),. = (%, 

We have evaluated including all terms of order (p1(pg, 
where a < 2 and 6 < 1. Therfore we can calculate the 
mean (q^) only up to order {p\ips- la principle with the 
Taylor series cut off in this way we can only evaluate 
Varga to this order as well. However, if one switches the 
variables to Qa = 53 — then it turns out that the Taylor 
expansion of Qa in and (ps contains no terms zeroeth 
order in tps ■ This fact allows us to compute Var Qa and 
hence Varga to order 

We now calculate (ga)^ to first order in (ps and pL- 


('?3)/i — 2 + -^0 


1 „_2(-/(| + i‘^L 


81 


vi)Vx<Ps)> 


{Jh 


(97) 


Now if we are dropping all terms second order in ps, 
then J can be considered a function of pL for the pur¬ 
poses of evaluating the numerator. Thus the numerator 
in Eq. llilTIl is the mean value of X'^ps times a function 
of pl- Since we are working only to order p\, and we 
know that {p^X'^ps)x = 0 , it follows that the numerator 
in Eq. (EU vanishes. Therefore 

{Q3)ti = ^ + 0{kL,pl,pl). (98) 


We will compute Var ga to first order in the power spec¬ 
trum of the long-wavelength modes or to second order in 
the long-wavelength modes themselves, i.e. to order p'j^. 
We will also work to second order in ps- Using Eq. P^. 
and defining Qa = ga — we find 

Varga = {QD^ - (Qa)^ = (Q^)^ _ 0{kl,p%,pl). (99) 

That is, the term (Qa)^ cannot contribute at order 
even though we have only computed (Qa)^ to order ps- 
Then 

^ {JQ^x ^•1nn^ 

Varga = —. (100) 


(<7)x — {^-\-^l{PL + Ps)~'rV2{p^‘L + ‘2-PLPs))x — + 

( 101 ) 

The numerator is 

/TZ. 2 X 25/3 , 625 / 3 . 2 /-I , 2500 / 3 . 1 /-I 

1 U'a/x + I45g7j4 + 729iJ() 

^ 5000/3/_i , 1251 / 1 / 3 ,i/_i 
729i/4 243i/4 

, 500:/i/3/-i , 2bv2hl-i 
+ 243//4 + 817/4 • 

The physical origin of the terms in Eq. nu^ii is as follows: 
the first term is the usual variance of the deceleration 
due to some patches of the universe being over- or under- 
dense. The second term is the “smearing” loop correction 
described above. The third term is a correlation between 
the order X'^ps and plX'^ps in ga fEo. ITTIl that arises 
from the nonvanishing bispectrum from inflation. The 
fourth term comes from two places: the variance of the 
PlX'^Ps term (which is the infrared-divergent term iden¬ 
tified by BMR) and the correlation between X'^ps and 
p\X'^ps terms in Eq. llTTIl . In this case, the variance 
term contributes a coefficient of and the correlation 
contributes yielding the total coefficient of 

The fifth, sixth, and seventh terms represent the modu¬ 
lation of the earlier terms by the non-coordinate measure. 

Looking at Eqs. (irmTi and (ITITa . it appears at first 
glance that /_i affects Varga in a highly nontrivial way 
depending on the values of /a, Zap, and 1^ 2 - This is in 
fact not the case, because /a. Zap, and /a ,2 are not inde¬ 
pendent quantities. Integration by parts in the definition 
(Eq.El and dropping boundary terms gives 

^,^ = -{■0+1)1^,(3-1. (103) 

Noting that the definitions imply /a = Za^o, we conclude 
that Zap = — 4/3 and hp = I 6 Z 3 . One might object that 
the boundary terms cannot be neglected because Z 3 is 
ultraviolet divergent. Of course, if this divergence is not 
regulated, Z 3 = 00 , in which case even the first-order per¬ 
turbation theory result for Var ga is infinite and it makes 
little sense to talk about higher-order corrections. Physi¬ 
cally, the divergence is regulated by putting in some cut¬ 
off in wavenumber k. In the case of KMNR, the regulator 
is placed at the horizon scale, k ~ Hq, where one sepa¬ 
rates super-Hubble perturbations from the sub-Hubble 
modes (the latter are considered to be observable per¬ 
turbations to the universe, rather than a correction to 
the observed scale factor). In the case of deceleration 
parameter measurements from SNe la, the cutoff scale 
is roughly of the order of the distance D to the super¬ 
novae; modes with wavelengths less than this cannot be 
treated by a local deceleration parameter based on the 
second derivative of the d^-z relation at z = 0. In either 
case, the cutoff occurs at a fixed physical scale (the hori¬ 
zon scale or the distance to the supernovae), rather than 
at a fixed coordinate scale x. Therefore the position of 
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the cutoff in terms of coordinate wavenumber k is mod¬ 
ulated by (fL as ^cutoff oc Since this is the same 

as the modulation of the wavenumber in the inflationary 
power spectrum (k) by large-scale modes (see Eq. , 
one can simply absorb the cutoff into P!f\k). Then 
becomes finite and Eq. is valid. 

Substituting our results that /a i = — 4/3 and I 3 2 = 
I 6/3 into Eq. (I102II reduces it to 


{JQ\)x — 


25/3 

81P4 


(1 -h V2I-1) 


and hence 


Var ga = 


25/3 

81P4’ 


(104) 


(105) 


independently of the value of /_i. The long-wavelength 
contribution to /_i, which is responsible for the infrared- 
divergent terms, cancels out when all contributions are 
included. This is true for both the coordinate and the 
particle number measure. This specific example confirms 
the general conclusions based on causality found above. 


D. Interpretation 

The correspondance between our results and BMR re¬ 
veals the origin of the infrared divergence: it comes from 
the assumption that (ps was taken as a function of x 
rather than y. Since the x coordinate system is dis¬ 
torted by superhorizon perturbations, these superhorizon 
modes distort structures within the horizon and affect the 
luminosity distance. Whether this effect is observable de¬ 
pends on whether the values of ps within the horizon are 
a Gaussian random field with the usual power spectrum 
Pip{k) in X or in y. This can be answered only by a theory 
of the initial conditions. If the initial conditions are set 
by single-field inflation, then the fluctuations ipL that de¬ 
termine the relation between the coordinates x and y are 
set down when these large scales leave the horizon. They 
then become classical, and later on (i.e. many e-folds of 
inflation later) the perturbations ips are generated. Since 
the generation of ps rnust be causal, one would expect 
that within regions small compared to the wavelength of 
Pl, inflation generates ps homogeneous and isotropic in 
the y coordinate system. If the calculation is done in 
the X coordinate system, as we did in Sec. ITv7!l then 
the infrared divergence from the second-order perturba¬ 
tion theory found by BMR cancels against three other 
infrared divergences: one arising from the correlation of 
first- and third-order perturbation terms, one from the 
correlation of first- and second-order perturbation terms 
that arises from the primordial bispectrum, and one from 
loop corrections to the power spectrum predicted by in¬ 
flation. Physically, the disappearance of /_i from the 
statistical properties of observables such as Var 53 is a 
manifestation of the fact that inflation wipes out initial 
conditions: the later stages of inflation prevent one from 


observing the pre-existing larger-scale structure of the 
universe, including the perturbations generated during 
the early stages of inflation. (Recall that this is also how 
inflation solves the flatness and homogeneity problems.) 

The scenario proposed by KMNR is closely related to 
the infrared divergence. Specifically, they argue that the 
gi 0 (pi:,/ 3 y 2 term in Eq. (I74II has a large variance that 
causes the deceleration parameter to also have a large 
variance. At second order in perturbation theory, this 
term is the infrared-divergent pL^/l-ps found by BMR, 
which we find to be canceled by other divergences that 
were not considered by BMR. KMNR uses the full pref¬ 
actor hence includes the divergence associ¬ 

ated with the p'j^V^ps term, but does not consider the 
divergences from the bispectrum or the loop correction 
to the power spectrum. Erom the perturbative calcula¬ 
tion of Sec. IWT!1 it is not obvious whether this cancel¬ 
lation extends to arbitrary order. However if the later 
e-folds of inflation produce perturbations homogeneous 
and isotropic with power spectrum P^\k) in the y co¬ 
ordinate system, as must happen on account of causality, 
then we need only consider Eq. m to realize that the 
super horizon structure is irrelevant, and the cancellation 
of the infrared divergences that arise in the x coordinate 
system must be exact. 

[As noted in Sec. El KMNR use qi rather than as 
their deceleration parameter. As far as Eq. m is con¬ 
cerned, this does not matter since Eq. IZD is valid for 
qi as well as for 93 . To see this, note that in first order 
perturbation theory, if one does a spherical expansion of 
the perturbation around the observer, symmetry implies 
that only the I = 0 modes contribute to qi or 93 . If 
only the I = 0 modes are present then we have 53 / in¬ 
dependent of direction n* and 93 = 94 . Since Hi = H/i, 
a comparison of Eqs. 03 and 1123 with the assumption 
that cr = w = 0 from spherical symmetry implies qi = 94 . 
We then have qi = 93 by transitivity. Hence Eq. 03 and 
thus Eq. iTTHi are valid for qi instead of 53 .] 


V. CONCLUSIONS 

In this paper, we have investigated the KMNR ex¬ 
planation for the accelerating expansion of the universe, 
which suggests that (i) in sufficiently underdense regions 
of the universe, the Hubble expansion appears to accel¬ 
erate (gi < 0 ) even with only normal matter present and 
Einstein gravity; and (ii) the variance Var qi of the decel¬ 
eration parameter is much greater than the simple cal¬ 
culation ~ ( 10 “^)^ because of the influence of perturba¬ 
tion modes with wavelengths many orders of magnitude 
larger than the Hubble length, so that the acceleration 
gi < 0 has a non-negligible probability of actually occur¬ 
ring. We have shown that suggestion (i) is ruled out by 
the Raychaudhuri equation, i.e. with GR and the SEC 
one always has qi > 0. We have also shown that sug¬ 
gestion (ii) is not true for the perturbations generated 
by single-field inflation, which does indeed predict Var qi 
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of order /| ~ (assuming an ultraviolet cutoff at 

the Hubble scale, as in KMNR). Therefore the KMNR 
mechanism cannot explain the acclerating expansion of 
the universe. 

What freedom is there to construct models related to 
KMNR that do explain the acceleration of the universe 
with superhorizon perturbations? Of the two elements of 
the KMNR mechanism, (ii) may be the easiest to circum¬ 
vent. In single-field inflation, the perturbations laid down 
by the early stages of inflation are adiabatic so the per¬ 
turbations laid down by the later stages of inflation must 
“look the same” (i.e. have the same power spectrum in 
our y coordinate system) everywhere. This requirement 
is the physics underlying the particular limiting forms 
for the bispectrum and the one-loop correction to the 
power spectrum that we used in Sec. ITW 1 to argue that 
single-field inflation produces no infrared divergences. In 
contrast, multi-field inflation models can produce isocur¬ 
vature modes so that different patches of the universe 
look different. In this case Vargi may be large because 
one takes the variance over regions with a different mix of 
cosmological fluids |^. Even then, however, an observer 
can only see one Hubble volume because of causality and 
so Var qi is not an observable, rather it is the variance 
of a distribution from which one obtains a single sample. 
In such a case the only role played by the superhorizon 
perturbations is to alter the initial conditions: all observ¬ 
ables, including 51 , are obtained from causal evolution 
that can be calculated from knowledge of the perturba¬ 
tions within the observer’s horizon. More importantly 
for the dark energy question, the combination of GR and 
the SEC still forbids gi < 0. 

If one is to find a way to keep GR, the SEC, zero vor- 
ticity, and neglect perturbations with wavelengths small 
compared to the Hubble length, and still maintain consis¬ 
tency with the observational results, one is forced to find 
a situation in which <71 as defined above is not what is ac¬ 
tually measured in the SN la experiments. BMR argued 
that <73 is a better representation of what is observed, 
since it is based on luminosity distances; the concordance 
cosmology has (73 ^ —0.6. Unfortunately, we found in 
Sec.^that (73 < 0 is possible within the GR-I-SEC frame¬ 
work only for anisotropic expansion, since the angular 
average (i? 3 /( 73 ') 47 r must be non-negative. In this case, 
there must also be lines of sight along which ( 73 / > 0 . 
There is no observational evidence that this is the case. 
Indeed, it would be an extraordinary mystery to have 
the deceleration parameter ( 73 / vary by an amount of or¬ 
der unity on scales comparable to the Hubble length and 
still produce a CMB isotropic to the level of a few parts 
in 10 ^ 

We conclude that cosmological models based on GR 
with irrotational initial conditions and perturbations 
only at and above the Hubble scale and only matter fields 
that conform to the SEG cannot explain the accelerating 
expansion. This paper does not exclude the possibility 
of using backreaction from sub-Hubble perturbations to 
explain the accelerating expansion. The latter possibility 


is difficult to investigate as it involves complicated non¬ 
linear physics; the perturbative calculations [^. which 
account for the nonlinear evolution of density perturba¬ 
tions, but not for strong field GR effects, suggest that 
the sub-Hubble backreaction is small |^. Nevertheless, 
only a full non-perturbative analysis would be definitive 
since there are rare regions of the universe such as black 
holes that cannot be described as a perturbation of an 
FRW spacetime. Regardless of the sub-Hubble physics, 
however, the superhorizon perturbations are not a viable 
mechanism to explain the acceleration of the universe: 
evolution of perturbations that lie beyond the horizon in 
real space cannot affect observables, and the superhori¬ 
zon perturbations can only act via their effects on the 
initial conditions within our horizon. These effects are 
then constrained by the “no go” theorems that require 
qi and (74 to be non-negative unless one invokes vorticity, 
modified gravity, or dark energy. 
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APPENDIX A: CONVERGENCE RADIUS OF 
SERIES FOR 4- 


The purpose of this appendix is to investigate the con¬ 
vergence properties of the power series Eq. for dt in 
the model of Sec. EH Our framework will be the result 
from complex analysis that a Taylor series of an analytic 
function has a radius of convergence given by the dis¬ 
tance to the nearest singularity. We begin by defining 
the dimensionless scale factor 


/3 


^0 


so that Eq. lUSl) becomes 



Oopen = ^(^^/^(cosh?? - 1) 


and 


= -(sinhr; — rf). 


(Al) 


(A2) 


(A3) 


It will be most convenient for the purposes of this ap¬ 
pendix to treat Eq. as a Taylor series in the dimen¬ 
sionless variable f3 rather than a or C. If we determine 
the radius of convergence in (3, then this immediately 
yields the radius of convergence in a or C. 

We want to determine the singularities of the func¬ 
tion dt of Eq. iBHIl or equivalently ln(aopen/a). First we 
comment on the analytical structure near /3 = 0. The 
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function ry(/ 3 ) defined by Eq. iTOll is double-valued near 
zero, with two solutions rji and —r]i. This presents no 
problem because both result in a single value for Oopen- 
Specihcally, near P = t] = 0 , we have |(sinh77 —77) ~ ^77^ 
so 77 ^ 2 / 3 ^/^ and 

Uopen ~ ^ 2EC^/^P - i^/3. (A4) 


Thus Oopen IS well-behaved at the origin (even for complex 
values). 

Equation IS 3 ) is the Taylor expansion of Eq. ll^ in 
P and hence the value of P with smallest absolute value 
that makes ln(aopen/a) singular determines the radius of 
convergence of Eq. lEUi. We can find these singularities 
by taking the derivative, 


d aopen(/3) 

dp a{P) 


d cosh 77 — 1 

sinh rj drj 
cosh 77 — 1 (i /3 


re¬ 


using the implicit derivative of Eq. foil . 


(A 5 ) 


only apparent since at small values of / 3 , we know that 
ln(aopen/a) 0 . Thus the only singularities can appear 
when cosh 77 = 1 and P ^ 0 . The solutions to cosh 77 = 1 
are 77 = — 2 'Kim where m is any integer. Any point with 
77 = — 2777777 , 777 7!^ 0 must actually be a singularity since 
for 77 = —2777777 -|- e, we have (by periodicity of the hyper¬ 
bolic functions with period 2777) 


2 sinh 77 


2 sinh e 


2e 


(cosh 77—1)2 (cosh £—1)2 (e^/2)2 


8e 


-3. 


(A8) 


then since / 3 ^/^ and /3 ^ are analytic for /3 0 there must 

be a singularity in Eq. The points with 77 = —2777777, 
777 0 correspond to 

/ 3 ^^^ = ^[sinh(—2777777) -I- 2777777] = ^ . (A 9 ) 

The values m = 0 correspond to /3 = 0 where there is no 
singularity. The other singularities occur at the points 


2/3 

/3 = I ^77 1 


(i-) 


(AlO) 





we find 


d , flopen _ 2 sinh 77 0I/2 0—1 

dp a (cosh 77 —1)2 


(A 6 ) 


(A 7 ) 


This function appears to have singularities where 
cosh 77 = 1 or /3 = 0 . The singularity at /3 = 0 is 


777 y^ 0 , where E is any cube root of — 1 . The closest such 
singularities to the origin have m = ±1 and hence |/ 3 | = 
(§77)^/^. Thus we conclude that the radius of convergence 
of Eq. is \P\ < (§77)^/^ or 

\ 2/3 

\CH^^a\ < i^-nj . (All) 
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